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^ , The phase structure of Nambu-Jona-Lasinio model with N-component 

. fermions in curved space-time is studied in the leading order of the 

1/N expansion. The effective potential for composite operator ipip is 
calculated by using the normal coordinate expansion in the Schwinger 
proper-time method. The existence of the first-order phase transition 
caused by the change of the space-time curvature is confirmed and the 
dynamical mass of the fermion is calculated as a simultaneous function 
of the curvature and the four-fermion coupling constant. The phase 
diagram in the curvature and the coupling constant is obtained. 
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It is quite interesting to make investigations of how the phase transition in quantum 
field theory takes place under the circumstance of the early universe. In paticular 
much interest has been taken in clarifying the mechanism of the spontaneous sym- 
metry breaking in the grand unified theories under the infiuence of the temperature, 
density and external gravity |]T|. 

In the standard grand unified theories the Higgs field appears as an elemen- 
tary field and has a Yukawa coupling to fermion fields giving masses to the fermions. 
As in the technicolor model for electroweak theory it is possible to assume that the 
Higgs field may be a composite system of some fundamental fermion fields. We shall 
call the model based on such an assumption generically the composite Higgs model. 

As far as we know a very little number of works have been reported in the 
study of the composite Higgs models under the circumstance of the early universe 0. 
Therefore we launched our plan to make a systematic study of the composite Higgs 
models under the circumstance of the early universe. Before going into a direct inves- 
tigation of the composite Higgs models it may be better to start with a fundamental 
study of prototype models of the composite Higgs particle with finite temperature, 
finite density and external gravity. In the present communication we take Nambu- 
Jona-Lasinio model as one of such prototype models and make a brief report of our 
work in the study of the phase structure of Nambu-Jona-Lasinio model in external 
gravity. 

The Nambu-Jona-Lasinio model in curved space-time is defined by the ac- 
tion, 



where g is the determinant of the space-time metric g^j_p, lti{x) the Dirac matrix in 
curved space-time, V^ip the covariant derivative of the fermion field ip, N the number 
of the fermion species. We work in the scheme of the 1/N expansion and perform 
our calculation in the leading order of the expansion. Our notation is basically in 



In practical calculations it is more convenient to introduce auxiliary fields a{x) and 





conformity with the ( 



) convention in the book by Misner, Thorne and Wheeler [Q. 
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7t{x) and start with the following action instead of Eq. (1 



S = d XV— ^ 



(2) 



In the flat space-time it is well-known that the global abelian chiral symmetry 
posessed by the Lagrangian (1) and (2) is broken spontaneously if the coupling con- 
stant A exceeds a critical value. Our purpose in the present paper is to see whether 
this phenomenon is modified under the influence of external gravity. We start with 
the generating functional given by 

Z[rj, v] = J DipDipDaD-K exp ^iS + if]ip + iiprj^ , (3) 

where 1] and fj are source functions. Integrating over the fermion fields with rj = f] = 
we find 

Z[0,0] = J DaDnexp{iNSeff} , (4) 
where we pulled out an obvious factor N in defining the semiclassical effective action 

^eff 

Seff = y'rf^Xv^|-^(a2 + 7r2)|-2lnDet07'^(x)V^-(a + 2757r)} . (5) 

In the leading order of the 1/N expansion the effective action F (with N factored out) 
is just equal to Sejf. Thus we obtain 

r[a,7r] = 5,^;[a,7r] + 0(l/N). (6) 

The effective potential (with N factored out) in the leading order of the 1/N expansion 
is then given by 

V{a, n) = ^(^2 + vr^) + iTT\n{x\{i-f>'{x)V ^ - (a + h,n)}\x) . (7) 

In Eq. (7) the variables a and vr are regarded as constant. 

To estimate the second term in the right-hand side of Eq. (7) we adopt the 
Schwinger proper time method^. For this purpose we rewrite Eq. (7) in the following 
form, 



V{(T, tt) = 77T-(o"^ + TT^) — iTi In S{x, x; s] 
2 a 



where the Green function defined by 



S{x,y;s) = {x\{iYV,-s)-'\ 



(9) 



is the solution of the equation 

{27''(x)V^ - s}S{x,y; s) 



-5\x-y). 



(10) 



We would like to calculate the Green function S{x,y,s) in the approximation of 
weakly varying gravity where we neglect any terms involving derivatives of the metric 
tensor higher than the third derivative. We use the Riemann normal coordinate 
expansion 0. Here we keep only terms independent of curvature R and terms linear 
in R. The calculation proceeds just in parallel with the one given in the paper by 
Parker and Toms[0. The result is 



S{x, X] s) 



(2^ 



{rka + S] 



1 



1 



+ -R^,k^k%rka + s. 



[k^-s^f 



where 



J 



ab 



(11) 



(12) 



and Latin indices with a caret symbol are vierbein indices. It may be seen from Eqs. 
(8) and (11) that V{a,7c) is symmetric in a and vr and so it is enough to discuss 
V{a, 0) insted of the full expression in the following arguments. Inserting Eq. (11) in 
Eq. (8) and making the Fourier transformation we obtain 

1 o 



Via,0) 



-a 



2A 
— zTr 



, [ d^k 
ds / - — -- 

J (27r)4 



i^k, + s 



1 



1 



+ -R^,k^k''{Yka + s)-^ 



k^ - s2 12 
1 1 



R{i% + s) 



1 



(A;2 - s2)2 



(P - s2)2 



(13) 
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To perform the momentum integration we make the Wick rotation and regularize the 
divergent integral by the cut-off method. We derive the following expression for the 
effective potential, 
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a'A^ + AMnl 1 + 



(47r2) 
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(47r)2 6 



a 
A2 



+ 



-aMn 1 + 



A^ 



AV^ 



A2 + (72 



(14) 



In Fig. 1 a typical behavior of the effective potential V{a, 0) as a function of 

curvature R is given for fixed four-fermion coupling constant A. As is well-known in 
the flat space-time, the chiral symmetry is broken if A > Aq with 

4^2 



An 



A2 



(15) 



In drawing Fig. 1 the coupling constant A is kept in the region A > Aq. It is clearly 
seen in Fig. 1 that the first-order phase transition takes place as R changes. Note, 
however, that the transition is found to be of second order for A < Aq . The dynamical 
mass of the fermion is calculated by analyzing the gap equation dV{a, 0) / da = : 

A2\ A2 



da 



a=ao 



'487r2 



R 



In 1 + 



+ 



+ 



A^ 



K^ + al (A2 + a2)2 







(16) 



The solution of the gap equation (16) corresponds to the vacuum expectation value 
of the composite field ipil) in the true vacuum and is equal to the dynamical mass of 
the fermion, m = ao. In Fig. 2 the dynamical mass of the fermion is plotted as 
a function of the curvature R for fixed A. The coupling constant A is kept again 
in the range A > Aq. The behavior of the dynamical mass as shown in Fig. 2 is 
characteristic for the relatively small couphng constant, Aq < A < 2Ao. For larger 
coupling, A > 2Ao, the behavior near the critical point R — Rcr is quite different 
from the one in Fig. 2: the curve representing the dynamical mass is bent upward 
near it! = Rcr- At any rate there is observed the gap in the dynamical mass at the 
critical point R — Rcr reflecting the nature of the flrst-order phase transition. By the 
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direct numerical analysis we find that there is no gap at i? = Rcr = if A < Aq and 
so the phase transition is of second order. 

It is possible to obtain the critical value of the curvature R and the coupling 
constant A by observing the behavior of the effective potential. The critical values Rcr 
and Acr constitute a critical curve in the R — X plane as shown in Fig. 3. It is found 
from Fig. 3 that for large positive R the chiral symmetry is restored even if A is kept 
in the region of the broken phase for R = 0. On the other hand the chiral symmetry 
is always broken for negative R irrespective of the value of A. This latter conclusion 
comes from the particular behavior of the effective potential (14) with vr = 0, i. e., 
the term proportional to R in Eq. (14) with vr = dominates over other terms if 
R < and makes the point cr = an unstable stationary point. The above peculiar 
behavior of our effective potential may be due to our approximation to keep only the 
first-order term in R in the normal coordinate expansion and should be left for the 
future detailed study. 

In the above discussions concerning Figs. 1, 2 and 3 we dealt with the rel- 
atively large value of the curvature, e. g., in Fig. 1 Rcr/A"^ = 0.656 which is large 
compared with ctq/A = 0.2 ~ 0.3 . Since we work in the normal coordinate expansion 
and keep only terms up to the one linear in R, large Rcr/A."^ compared with cr/A 
may spoil the validity of the expansion. One should, however, note that the terms 
quadratic or higher in R are not divergent by dimensional reason and the coefficient 
of R^ does not carry logarithmic functions of cr^/A^. Hence the R^ term is relatively 
suppressed compared with the term linear in R. Accordingly relatively large Rcr/ 
may be within the scope of our approximation. 

In the present model we found within our approximation that the first-order 
phase transition takes place as the curvature varies and the fermion mass is generated 
when the curvature is smaller than the critical value. Our conclusion has some simi- 
larity with the previous results in the massless scalar theory with cj)'^ couplingP|. The 
four- fermion theory in 2 dimensions (Gross-Neveu model ||^) has been investigated 
in curved space-time [^] and there the second-order phase transition is observed in 
contrast with our result. 
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FIGURE CAPTIONS 



Fig. 1. The typical behavior of the effective potential V is shown for fixed A 
(= 1.25Ao, Ao = 47r^/A^) as a function of the curvature it! where 
Rcr = 0.656A2. 

Fig. 2. The solution ctq of the gap equation (17) is shown as a function of the 
curvature R where A = 1.25Ao and Rcr = 0.656A^. 

Fig. 3. The phase diagram in A and R where Xcr and Rcr represent the critical 
values of A and R with which the region preserving the chiral symmetry 
is divided from the one of the broken chiral symmetry. 
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